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We explore QCD calculations for the process jp —>■ VX where V is a vector 
t^- \ meson, in the region s 3> —t and —t 3> Aq CD . We compare our calculations 

for the J/tp, <j) and p mesons with data from the ZEUS Collaboration at HERA 
and demonstrate that the BFKL approach is consistent with the data even 
for light mesons, whereas the two-gluon exchange approach is inadequate. 
We also predict the differential cross-sections for the T and ui for which no 



Oh 



X 



data are currently available. 

1 Introduction 



We study the process jp —> VX (where V is a vector meson: V = p, uj, 4>, J /ip or T) in the 
perturbative Regge limit, s>-t> Aq CD . The largeness of —t allows the application of 
the solutions of the non-forward BFKL equation. We apply the analytic solutions to the 
cross-section arrived at in jl|, |2| in the case of a delta-function distribution for the meson 
wavefunction. Specifically we apply this solution to the case of a real photon which has 
been measured at HERA [Q, |||. 

The cross-section for the helicity-flip process 7p — * VlX vanishes for a delta-function 
wavefunction and throughout this paper the cross-sections pertain to the process ^p — > 
VtX (where L and T correspond to longitudinal and transverse polarization respectively). 
This approximation is justifiable since the measured rate for longitudinal mesons is indeed 
small ||]. We make predictions and compare to experiment for the mesons for which data 
are available, the heavy J/ip and the light p and (f> (for which we might expect the delta- 
function to be a poor approximation). We show that the data for all the mesons can be 



well described by the BFKL approach. We then demonstrate that the two-gluon exchange 
model is incapable of adequately describing the data. Finally we make predictions for 
the T and to. 

The high momentum exchange allows us to factorize the cross-section into the usual 
product of the parton distribution functions and the parton level cross-section: 



^ = (&/(*>*) +?/(M)]J Jt +G(x,t) Jt , (1) 

where G(x,t) and qf(x,t) are the gluon and quark parton distribution functions respec- 
tively and we sum over flavour, /. For a large separation in rapidity between the parton 
and the vector meson we may write, 



da(jp^VX) ( 4iV c 4 \da(jq^Vq) 

where we define the parton level amplitude by 

^ ^IMf, (3) 

where s is the centre-of-mass energy squared of the photon-proton system and s' = xs is 
that of the photon-parton system. 



2 The Two Gluon and BFKL Amplitudes 

We consider two types of colour singlet exchange: two-gluon and BFKL. The two-gluon 
amplitude can be expressed to leading order in s in the impact factor representation: 

^■"° tf /p iCT - (4) 

where X 7 y and X qq are the impact factors associated with the couplings of the two gluons 
to the external particles. Vectors with the subscript _L are two dimensional transverse 
momenta and Q\ = —t. With the above definitions the impact factor for q —> q is given 
by 



Iqq = ~N c as ^ 

The 7 — > Vt impact factor, assuming a delta-function form for the meson wavefunction, 

is 



Cas 2N c [ Q\ + M 2 (k ± -Q ± /2) 2 + M 2 /4)> (6) 



where, 



C 2 = 3lf +e _ M v /a em . (7) 

My is the mass of the vector meson and Y\ e is the electronic decay width of the me- 
son. Using these impact factors we obtain the two-gluon amplitude as a function of one 
dimensionless parameter, r = —t/M v : 



qq \l-T 2 V At 
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which has the mathematical property J- Bom ->0asr->l. 

The BFKL amplitude, in the leading logarithm approximation (LLA), is given by ||, 

FbfklVJ) = -^jdu {u2 l\ /A)2 e^lf (Q ± )lf(Q±), (9) 

where 

X (v) = 4Ke(iP(l)-^(^ + iv\\ (10) 

and 

-=£K£> (11) 

In LLA, A is arbitrary (it need only be small compared to yV) and a s is a constant. 
The impact factors are used, in conjunction with the prescription of ||, to give 



d 2 k 
(2n) 



4 A (Q±) = J^lA(k ± ,Q ± )Jd 2 Pl d 2 P2 (12) 



P1P2 / \pV \pI 



Jkj_-pi+i(Q±-k±)-p2 



In the case of coupling to a colourless state only the first term in the square bracket 
survives since lA(k±,Q± = k±) = Xa(A;_i_ = 0, Q±) = in this case. After some work, 
one obtains 0: 



i q AQ±) = -^ qq 2- 2w \Q±\- 1+2lu l { t l ■ lv \ (13) 

r v, n v _ rr 7r 2 r(l/2-^) /Qjy- rV^Hoo du , Q± y +2u 

« W±J - L ^WQ3_r(l/2 + ii/)U ) h/2-ioo 2ni\M v /2) (Mj 

r(i - u- »i/)r(i - it + ^)r 2 (i/2 + u) 

r(l/2 + u/2 - ii//2)r(l - u/2 - ii//2)r(l - it/2 + iv/2)T(l/2 + u/2 + ii//2) ' 
Putting these into (M) we obtain, 



T 2 (l/2 + u)r(l - u - iu)T(l -u + iv) 
r(l/2 + -u/2 - zi//2)r(l - u/2 - iv/2)T{l - u/2 + iv/2)T(l/2 + u/2 + iv/2) ' 

3 Results and Discussion 

Using the full numerical calculation for the amplitude we convoluted the partonic cross- 
section for the mesons with the parton density functions of the protonF], integrating over 
x in the region 0.01 < x < 1 (for which the ZEUS data are quoted). We obtained fits 
to the ZEUS data in terms of three free parameters ||0. One parameter is a s and the 
others appear in the denominator of the logarithm defining the energy variable z, i.e. we 
take A 2 = (3M^ + j\t\. 

Initially we tried the simple prescription of |l|, ^] where j3 = 7 = 1. The J ftp has been 
fitted with a s = 0.20 for this parameterization previously JTJf] and as expected this gave 
a x 2 /dof < 1, however these parameters gave very poor fits to the p and <p. When we 
varied all three parameters independently we found that optimum fits took small values 
of 7 and in fact we were able to put this parameter to zero. 

Fig. [l] shows contours of constant x 2 . The contours continue beyond the region 0.1 < 
(3 < 10 without closing. Further imposing the requirement that s' > 10 j3My as a sensible 
criterion for the dominance of leading logs, a simultaneous fit for the three mesons rules 
out the region of a s > 0.2 in Fig. |]. In addition, sub-asymptotic leading log contributions 
become increasingly important as we increase the strong coupling to a s > 0.20 Q. Fig. | 
suggests that we can get acceptable fits for a < 0.17; while this is true we are then forced 

^rv-os LO §. 

2 Note that the differential cross-sections for the three mesons have an overall normalization uncertainty 
of approximately 10% (which cancels in the ratios of cross-sections). This normalization uncertainty 
does not significantly affect our results. 

3 H1 has recently performed a similar fit where they found a s — 0.22 |4|. 



to accept unnaturally small values of (3. We note that this range of 0.17 < a s < 0.20 is 
consistent with the value extracted from the Tevatron "gaps between jets" data [ 10 1 . 



It should be observed that the acceptable fits lie in a narrow valley in parameter space. 
It turns out that the relation, 

«■<» = raw? (16) 

with d = 0.358 gives a good approximation of the flow of a s with (3 for a constant \ 2 ■ 
The width of the x 2 /dof< 1 valley corresonds to a band 0.197 < a c < 0.204, with the 
minimum at approximately a c = 0.200. 

Perhaps the most natural fit is (a s , f3, 7) = (0.2, 1.0, 0.0). The BFKL differential cross- 
section predictions for the three mesons are shown in Fig. ^, along with the ZEUS data 

Our results suggest that the largeness of —t does allow a perturbative calculation to 
be performed however it is reasonable to ask whether it is necessary to invoke the full 
machinery of BFKL. Hence we also show the two-gluon exchange fits in Fig. ^. The two- 
gluon calculation only has one free parameter, a s , and the differential cross-section has 
a fourth order power dependence on it. It is apparent that this approximation provides 
a very poor description of the light mesons, for which r ~ 1 over much of the range of 
data. Note that the model starts to give the correct qualitative behaviour away from 
r = 1. Note in particular that a s has been fitted for each meson separately and also the 
unnatural trend that a s falls as My falls, which means that running the coupling only 
makes the situation worse. The J/ip differential cross-section, to which we would expect 
our calculation to be most applicable can be fitted by the two-gluon exchange amplitude 
in the ZEUS i-range. For a plausible strong coupling a s = 0.4 we obtain Xj/^/dof= 0.8. 
It should be noted however that the two-gluon curve is already starting its dip towards 
zero at \t\ = Mj,,, by the end of the ZEUS t range. HI data for this process reaching out 
to much higher \t\ should make it clear whether or not two-gluon exchange is sufficient 
to describe the J/ip- 

Fig. |H shows our predictions for the T and ui mesons. The bounds come from the 
uncertainty in the normalization of the LO BFKL solution, specifically from varying a s 
while keeping 0.1 < /? < 1.0 and x 2 /dof< 1.0. The major part of the uncertainty is in 
normalisation; the shape being rather well predicted. The bounds for the u are clearly 
much narrower than for the T. This is due to the fact that M p ~ M u which means that 
da^/dt ~ (re e /r^ e ) • d(Tp/dt. This scaling relationship between the p and w is a model 
independent feature and also true of the two-gluon model. It should be noted that as we 
approach the lower bound for the T in Fig. || (a), we are extending the model beyond 
the leading log region (due to largeness of /3M 2 -), so this bound must be viewed with 
some caution. 

If subsequent measurements yield data lying outside these bounds it would indicate either 
a failure of the LLA BFKL formalism, or a need to refine our treatment of the meson 



wavefunction. If other light vector mesons, such as the u), subsequently confirm our 
predictions it would suggest that in some instances even mesons made of light quarks 
can act as if they consist of two constituent quarks which share the meson's energy and 
momentum. 

For all of the measured mesons, the experimental errors are about 10-20% of the ab- 
solute value of the differential cross-section. That we can fit the data suggests that 
genuine higher order corrections and the corrections due to the meson wavefunction may 
contribute no more than 10-20%. 



Final Remarks 

In (ll]] a LLA BFKL calculation incorporating a more sophisticated meson wavefunction 
than that used here is performed. It was concluded in @ that this light meson calculation 
is incompatible with the available data. However, the author of |TT] interprets the masses 
appearing in (g) as current masses (as ought to be appropriate for a truly perturbative 
calculation) which justifies the neglect of the quark mass for light mesons, rendering |t| the 
only relevant scale. By interpreting the quark mass as a constituent mass, and assuming 
that the quark and anti-quark share the meson momentum, we have demonstrated that 
the data can be understood provided the constituent mass is taken to be the scale in 



A. Another point of deviation is in the treatment of the strong coupling a s . In |TT|] the 
coupling is a running coupling whilst we have considered a fixed coupling. The problem 
with running the coupling it is that we do not really know how it runs. Of course a 
complete description should incorporate a running coupling, but the work of [12| might 
be a hint that a fixed coupling may be appropriate for the BFKL exchange. The work 



of [11] has been developed in [|T^] where it is proposed that a large contribution arising 
from qq fluctuations in a chiral-odd spin configuration should play an important role in 
the region of the data. It remains to be seen if the simple model presented here can be 
justified within this approach. 

Our results suggest that the two-gluon model is inadequate at least for the light mesons. 
The model predicts a dip at r = 1 which is not present in the data and this is the biggest 
obstacle to getting a decent fit. Corrections to the two-gluon approximation do fill in 
the dip (see [[jj for more details on how this occurs) and further study is called for to 
establish whether a more sophisticated two-gluon or finite order (in a s ) calculation might 
work. 

The model explored here has experimental tests to face in the future. If data confirm 
our predictions for the T and u> it will be impressive. Data on the process *yp — ► yX, for 



which theoretical calculations have already been performed |14[| , should be obtained in 
the forseeable future. This process along with the ones considered here will continue to 
provide an important test of the validity of BFKL dynamics. 
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Figure 1: Contour plot of x 2 /dof where \ 2 is that of the J/ip, <fi, p mesons com- 
bined. From the inside heading outwards the contours are those of x 2 /dof = 
1, 3, 9, 27, with the shaded region corresponding to x 2 /dof< 1. 
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Figure 2: Comparison of BFKL (solid line) and two-gluon (dotted line) calculations com- 
pared to ZEUS data [[| for (a) J/vt, (b) cf> and (c) p meson production 0. The 
BFKL curves correspond to the fit (a s , (3, 7) = (0.20, 1.0, 0.0), described in 
the text. The two-gluon curves are obtained by optimising a s in each case, i.e. 
(a) a s = 0.40, (b) a s = 0.35 , (c) a s = 0.27. 
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Figure 3: Predictions for (a) T and (b) u meson production. The bounds come from 
varying a s with 0.1 < j3 < 1.0 and keeping x 2 /dof< 1. The upper and lower 
bounds correspond to (a s , (3) = (0.175, 0.1) and (a s , (3) = (0.197, 1.0) respec- 
tively. 
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